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[57] ABSTRACT 

A system for use in detecting earth crustal deformation 
using an RF interferometer technique for such purposes 
as earthquake predictive research and eventual opera- 
tional predictions. A lunar based RF transmission or 
transmissions from earth orbiting satellites are received 
at two locations on Earth, and a precise time dependent 
phase measurement is made of the RF signal as received 
at the two locations to determine two or three spatial 
parameters of the antenna relative positions. The re- 
ceived data are precisely time tagged and land-line 
routed to a central station for real-time phase compari- 
son and analysis. By monitoring the antenna relative 
positions over an extended period of months or years, 
crustal deformation of the earth can be detected. 

6 Claims, 2 Drawing Figures 
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SYSTEM FOR NEAR REAL-TIME CRUSTAL 
DEFORMATION MONITORING 

ORIGIN OF INVENTION 5 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-568 
(72 Stat. 435; 42 USC 2457). 10 

BACKGROUND OF THE INVENTION 

This invention relates to a method for near real-time 
crustal deformation monitoring using long-baseline 
radio interferometry. 15 

There has been a continued and growing interest 
using long baseline radio interferometry for monitoring 
crustal deformation of the earth for prediction of earth- 
quakes. See “Geophysical Applications of Long Base- 
line Radio Interferometry,” by I. L Shapiro and C. A. 20 
Knight in Earthquake Displacement Fields and the Rota- 
tion of the Earth, edited by L. Mansinha, et. al., and 
published by D. Reidel, Dordrecht, Netherlands, 1970. 

The present inventor has made a study of the feasibility 
of radio interferometry in this respect and the results 25 
were reported in a paper titled “Radio Interferometry 
for International Study of the Earthquake Mechanism”, 
Acta Astronantica, VOL. 1, pp. 1427-1444, Perpamon 
Press (1974). The technique is not unlike that employed 
in stellar interferometers. It is operable on a very wide 30 
range of baseline lengths from 100 m to intercontinental 
separations. The term “Astronomical Radio Interfero- 
metric Earth Surveying” has been adapted for the tech- 
nique, now often referred to by the acronym ARIES 
described in that report. 35 

A precision geodetic measurement system based on 
ARIES has been designed and implemented through 
the use of a 9-m transportable antenna and fixed anten- 
nas, such as the NASA 64-m antenna of the Deep Space 
Communication Complex at Goldstone, California, and 40 
the 40-m telescope of the Owens Valley Radio Observa- 
tory. The system requires extragalactic random (quasar) 
radio signals received by the two antennas on a very 
long baseline to be recorded on magnetic tape for off- 
line cross correlation. The system will detect subtle 45 
motions of the Earth’s crust in three dimensions. Such 
motions are believed to precede an earthquake. If that is 
so, then a reliable earthquake prediction technique 
could be developed. The key to the system is the capa- 
bility of measuring the difference in arrival of identical 50 
quasar radio signals at the two antennas on the baseline 
using atomic clocks as timing references and high speed 
computer techniques for cross correlation. The problem 
in adapting the system to real-time monitoring, or to 
near real-time monitoring, using quasar interferometry 55 
is cumbersome for real-time observations because of the 
need for utilizing wideband data communications lines 
and a plurality of quasar sources to eliminate the effects 
of system drifts and achieve full three dimensional mea- 
surements. 60 

SUMMARY OF THE INVENTION 

In accordance with the present invention, real-time 
crustal deformation monitoring can be achieved by 
independent stations in a radio interferometry system 65 
with a pair of antennas on earth illuminated by an arti- 
ficial narrow spectrum radio signal from a source on the 
moon. Telephone circuits are used for real-time trans- 
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mission of the radio signal received by the two antennas 
to a central station for cross correlation together with 
data from synchronized clocks at the stations for use in 
compensating for transmission delay differences in the 
telephone lines. 

BRIEF DESCRIPTION OF THE DRAWINGS 

FIG. 1 illustrates the conceptual elements for moni- 
toring stations on a very long baseline. 

FIG. 2 illustrates a spherical triangle which relates 
the right ascension, a, and declination, 6, of the source 
to the right ascension, a& and declination, 6^, of the 
direction of the baseline vector. 

DESCRIPTION OF PREFERRED EMBODIMENT 

Referring to the drawings, real-time crustal deforma- 
tion monitoring can be achieved by independent sta- 
tions in a radio interferometry technique with a pair of 
small (2.5m) antennas 10 and 11 illuminated by S-band 
signals, or other artificial narrow band radio signals 
from a source 12 on the moon. The Apollo Lunar Sur- 
face Experiments Package (ALSEP) units placed on the 
moon by Apollo Flights 12 , 14 , 15 , 16 and 17 represent 
a valuable source of radio illuminations for very long 
baseline interferometers (VLBI) on earth. Such interfer- 
ometry has been reported by several investigators. 
Counselman, C. C., III, Hinteregger, H. F., and Sha- 
piro, I. I., “Astronomical Applications of Differential 
Interferometry,” Science, Vol. 178 , pp. 607-608, 1972; 
Counselman, C. C., III, et al., “Lunar Baselines and 
Libration from Differential VLBI Observations of AL- 
SEPS ” The Moon, Vol. 8, pp. 484-489, 1973; and 
Slade, M. A., MacDoran, P. F., and Thomas, J. B., 
“Very Long Baseline Interferometry (VLBI) Possibili- 
ties for Lunar Study,” in The Deep Space Network, 
Technical Report 32-1526, Vol. XII, pp. 35-39, Jet 
Propulsion Laboratory, Pasadena, Calif., 1972. How- 
ever, no use has been made of the ALSEP signals for 
terrestrial geodesy, much less for crustal deformation 
monitoring in near real time. 

Although the ALSEP signals will allow only fringe 
frequency observations and measurement of equatorial 
baseline components, 82% of vertical deformation can 
be sensed at a 34° latitude. This system will occupy sites 
previously measured in three dimensions by ARIES and 
monitored thereafter by this small antenna system. By 
comparing the ARIES equatorial baseline components, 
derived from quasar positions, with those determined 
by ALSEP signals, it will also be possible to relate lunar 
motions relative to the quasar reference frame within 
approximately 0. 1 arc sec. 

Unlike the quasars which emit very wideband radio 
spectral noise, the ALSEP is spectrally narrow, it being 
composed of telemetry signals from lunar-based scien- 
tific instruments. The ALSEP units are not radio tran- 
sponders of the type used aboard interplanetary space- 
craft. The ALSEP units do accept uplink commands 
from earth to control certain functions. However, no 
coherent response occurs from the device. 

Although the ALSEP spectrum is narrow, about 6 
KHz at 3 db, the radio flux arriving at the earth is equiv- 
alent to 50,000 Jansky (10“ 26 w/m 2 /Hz). By compari- 
son, in ARIES experiments a strong quasar is consid- 
ered to be a source of 5 or more Jy. This factor of 10 4 
difference in flux is then available to be traded off to 
achieve a simpler radio interferometry system, but at 
the expense of having introduced radio source dynam- 
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ics into the problem. Fortunately, LURE (Lunar Rang- 
ing Experiment) and the ALSEP/Quasar VLB! experi- 
ments, have been successfully dealing with such prob- 
lems and now represent a valuable resource for achiev- 
ing the desired radio interferometric geodesy simplifica- 5 
tions. See Bender, P. L., et ah, “The Lunar Laser Rang- 
ing Experiment,” Science, Vo! 182, pp. 229-238, 1973, 
and Slade, M. A., Preston, R. A., Harris, A. W., 
Skjerve, L. J. f and Spitzmesser, D. J., “AL3EP - Quasar 
Differential VLBI”, in the Deep Space Network jq 
P rogress Report, 42-33, Vol 2, pp. 37-47, Jet Propul- 
sion Laboratory, Pasadena, Calif., 1976. 

The RF spectral narrowness of the ALSEP radio 
signal does represent a limitation relative to the ex- 
tremely wide RF continuum emissions of quasars. A ^ 
broad RF emission spectrum allows the interferometer 
to precisely measure the delay function by signal cross- 
correlation within 0. 1 nanosecond, but a narrow spec- 
trum allows only the time rate of delay change to be 
precisely measured. See Williams J. G., “Very Long 2 q 
B aseline Interferometry and Its Sensitivity to Geophys- 
ical and Astronomical Effects,” The Deep Space Net- 
work, Space Programs Summary, Tech Rep. 37-62, 
Vol. 2, p. 49, Jet Propulsion Laboratory, Pasadena, 
Calif., 1970. Three dimensional geodesy requires inter- ^ 
ferometry delay measurements from at least two qua- 
sars at substantially different declinations. In actual 
practice, ARIES experiments use ten to twenty quasars 
observed over periods of 8 to 26 hours because of a need 
for an alternate solution strategy caused by frequency 
system instabilities. 

Having access to only delay rate data, also called 
fringe frequency, results in the interferometer being 
sensitive to only the two equatorial components (X and 
Y) of the baseline vector. The Z component insensitiv- ^ 
ity results in a loss of 18% of the information and some 
uniqueness of interpretation. More specifically, con- 
sider the effect of applying the rotation matrix (valid for 
latitude 35° ) to transform a unit local vertical displace- 
ment into geocentric baseline components: 
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Thus, a one meter increment in the local vertical ap- 
pears as a 39 cm decrease in the geocentric X compo- 
nent (direction toward Greenwich) and a 73 cm de- 
crease in the Y component (90° E of Greenwich). The 
spin axis component, which will not be sensed by this 50 
system, experiences a 56 cm increase. 

The narrow RF spectrum of ALSEP does offer an 
important advantage of being easily recorded or trans- 
ferred via telephone lines 14 and 16, A 3 KHz informa- 
tion bandwidth is suitable from signal to noise ratio 55 
considerations and quite reasonable for telephone cir- 
cuit transmission for real-time cross-correlation signal 
detection processing. 

The phone lines have instabilities that need frequent 
calibration. This can be accomplished by timeformat- 60 
ting the data in the usual Very Long Baseline Interfer- 
ometry (VLBI) manner at each station using clocks 18 
and 20. Each of the antenna stations will have a cesium 
clock which can be synchronized to the other clock 
within one microsecond or better, the time at each sta- 65 
tion is continually transmitted as spurs (pulses) over the 
telephone lines with the ALSEP signal, and the time 
received from each station is compared at a cross-cor- 
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relator station 22. Thus, the phone line delay can be 
determined well within the 3 KHz bandwidth restric- 
tions of the line. 

Consider two or more 2.5 meter (~ 8 ft.) diameter 
dish antennas with efficiencies of 50%, 200 Kelvin S- 
band receiver operating temperatures, 3 KHz band- 
width, 10 second coherent integration time and receiv- 
ing ALSEP signals equivalent to 5X 10 4 Jy. The signal 
to noise ratio (SNR) of the cross-correlation output is 
then 27. See the paper by the present inventor in Acta 
Astro nautica, supra. 

By deriving the interferometer phase every 10 sec- 
onds over an 800 second interval it will be possible to 
extract the fringe frequency with a precision of 50 to 
100 /xKz of S-band. Several 800 sec. samples of the 
fringe frequency over a wide range of sky visibility for 
the ALSEP units will allow the amplitude and phase of 
the diurnal signature to be determined along with the 
frequency system offset. Thus, equatorial baseline com- 
ponents can be measured. See Williams, supra, and the 
paper by the present inventor, supra. Having 50 /xHz 
data quality every 800 sec. implies 10 cm baseline preci- 
sion. 

The independently operated atomic frequency sys- 
tems at each station introduce another error source. 
Ideally, a frequency system such as a hydrogen maser 
could be used with virtually no contribution to the 
allowable error. However, cost aspects make it impor- 
tant to consider other frequency sources for this study 
where unattended remote station operations are highly 
desirable. At 800 sec. a cesium frequency system having 
a stability of about A f/f= 1.6 X 10 - 13 at S-band causes a 
350 fxHz noise level. Incoherently combining these 350 
jitHz observations over a 10 hour observation period 
allows improvement to an effective 50 /xHz observation 
which is equivalent to about a 10 cm baseline precision. 

Because the equatorial baseline information is inher- 
ently diurnal in nature, special care must be taken to 
either stabilize the receiving system or to explicitly 
measure possible thermal phase variations over the 
tracking pass. Such phase calibration can be performed 
and transmitted along with the time-formatted video 
data. 

Based upon LURE data presented by Williams, J. G., 
“Lunar Laser Ranging: Present Results and Future 
Plans,” EOS, Trans., Am. Geophys. Union, Vol. 56, 
970, 1975, the lunar ephermerides are believed to be 
internally consistent to within approximately 0.01 arc 
sec. Considering the two small antenna stations to be 
deployed 100 Km apart at established ARIES sites, a 
0.01 arc sec lunar position internal error causes a 0.5 cm 
baseline error. Possible angular rotations between the 
lunar frame and the quasar frame might be present and 
could be sensed at about the 0. 1 arc sec level by this 100 
Km baseline. A 100 Kin baseline is short enough that 
diurnal variation of the ionosphere are seif-cancelling in 
the differential measurement that is intrinsic to interfer- 
ometry. Residual effects due to ionospheric clouds of 
scale < 100 Km are estimated to be 10 cm per pass and 
random on a pass to pass basis. 

It is also possible to transmit surface meteorology 
through the phone line. The radiosonde calibrations 
will be accurate to 3 cm at zenith and considering that 
tracking data will be taken to only 20° elevation angles, 
the baseline will be affected at about the 10 cm level on 
a single pass basis. Water vapor radiometers could cali- 
brate the set tropospheric delay effects with a path 


4,170,776 


5 

equivalent accuracy of 2 cm, however, such instru- 
ments are relatively expensive and may be non-optimal 
for continuous monitoring on a time scale of months to 
years. Solar hydrometers may serve the wet tropo- 
spheric calibration need. 

In summary, a radio interferometry method has been 
proposed which uses the relatively strong S-band radio 
signals coming from the ALSEP stations on the moon. 
These strong signals allow several simplifications to 
independent station radio interferometry such as small 
2.5 m antenna sizes and real time cross correlation. The 
frequency system requirements of Af/f~2x 10- 13 , 
r=800 sec, Are obtainable commercially. Stations of 
the type discussed could be implemented and operated 
at comparatively low cost and could make important 
geophysical contributions to the determination of 
crustal deformations. Commercial telephone lines cou- 
pling the antennas stations to the cross correlator make 
near real-time cross correlation possible. The cross-cor- 
relation data are then used with calibration data 24 by a 
computer 26 programmed for equatorial baseline com- 
ponent estimation and crustal deformation monitoring. 

The sensitivity to geophysical and astronomical ef- 
fects have been investigated by J. G. Williams, JPL 
Space Programs Summary 37-62, Vol II, pp. 49-55, 
Mar. 31, 1970. Although the system then contemplated 
necessarily involved recording the output of two an- 
tenna separately on magnetic tapes later brought to- 
gether for processing, the geophysical and astronomical 
effects remain essentially the same where, as in the 30 
present invention, the output of the two antenna are 
coupled to a cross correlator for near real-time process- 
ing by conventional telephone. By time tagging the 
separate outputs with synchronized clock outputs, it is 
possible to tolerate some real time delay in transmission 35 
from one antenna longer than from the other, and to 
even experience changes in transmission delays as the 
telephone network rerouts the output of one antenna 
according to requirements of ever changing traffic on 
the telephone lines. By time formatting the separate 40 
data streams with the synchronized clock outputs, it is 
possible to virtually eliminate the delay effects of the 
telephone transmission lines linking the receiving sta- 
tions to the cross correlator. Consider a 6 KB/s data 
transmission rate corresponding to a 3kHz information 45 
bandwidth. A frame of data might consist of 600 bits 
with the first 6 bits all binary l’s and the first of the 6 bits 
synchronized to a 10 pulse per second output from the 
clock. The remaining 594 bits of each data frame would 
consist of a binary sampling of the station receiver but- 50 
put which would appear as an almost random series of 
binary bits. In this manner it would be possible to con- 
struct at the cross correlator a software algorithm 
which synchronizes the two arriving data streams every 
0.1 seconds without regard to telephone variatons up to 55 
approximately 50 milliseconds. To guard against delay 
and variatons greater than 50 m sec, a 10 bit series of all 
l’s might be used as a beginning of frame code every 
integer second with leading bit coincident with the one 
pulse per second output of the synchronized clock. 

At the cross correlator, the arriving bit stream from 
the station geometrically nearest to the moon will be 
buffer delayed typically by 5 bits or less corresponding 
to 0.84 m sec (allowing station separations up to 250 
Km). Adjusting the bit streams plus or minus 5 bits will 65 
keep the cross correlation function maximized so that 
the time variable nature of the cross correlation func- 
tion (i.e. fringe frequency) can be sensed and analyzed 
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for the interferometer equatorial baseline information 
being sought. 

The technique of interferometry as investigated by 
Williams, supra, but applied to the present invention, 
5 will now be briefly reviewed. Assuming the two anten- 
nas are separated by a distance D. The direction of the 
radio source 12 makes an angle i// to the baseline be- 
tween the antennas. There will be a difference with time 
of reception of the signal at the two antennas propor- 
10 tional to the angle i/j, at the farthest antenna from the 
source. The time delay r gf is given by 

T g =D/Ccos if/ (1) 

15 where c is the speed of light and D is the baseline dis- 
tance. The electric field at antenna 11 is given by 

E\ i —A cos a it (2) 

20 where A is the amplitude, o> the frequency, and t the 
time. The electric field of the other antenna is given by 

E\o=A cos a) (l— T g ) (3) 

25 The voltages received by the cross correlator are given 
by 

V\\{t)=A cos o/t+fyll (4) 

^10(0 = ^ cos ttf— CdTg+<|>10 (5) 

Phase shifts <J>11 and (f>10 are introduced by instrumen- 
tal and atmospheric delays which are explicitly cali- 
brated (24) in the data analysis. 

The frequency of the signal from the antenna 10 is 

F=o) — o)dTg/dt (6) 

The quantity of ct>drg/dt is known as the fringe rate vf- 
It is the negative doppler shift between the two an- 
tenna sites. Cross correlaton allows the time delay r g 
and the fringe rate vf to be measured. 

To be useful, the angle if/ must be related to the right 
ascension and declination of the source. FIG. 2 illus- 
trates a spherical triangle which relates the right ascen- 
sion, aand declination, 8, to the right ascension, at, and 
declination, 8b, of the direction of the baseline vector. 
From spherical trigonometry, the law of cosines gives 

cos ij/=sin 8 sin E/,+cos 6 cos S^cos (a*— a) (6) 

The baseline is fixed to the earth so that its right 
ascension ab can be related to its constant longitude \b 
and the right ascension of Greenwich a g (t) by 

aA=X 6 +a g (r) (7) 

Substituting Eq. (6) in Eq. (1) gives 

T = Tg+T 0 

60 where r 0 represents the synchronization of the station 
clocks, 

Tg~D/c[ sin 8 sin 8^+cos 8 cos S^cos (a*— a] (8) 

The fringe rate is given by 

vp— — aD / c oi e cos 8 cos db sin (a a)+v 0 (9) 


where o) e is the rotation rate of the earth, 0.73 x 10~ 4 
rad/s and v 0 is the frequency offset caused by the fact 
that the clocks run at slightly different rates. If v/? is 
expressed in Hertz, then 


v/r=— Z?w t </ Acos 6 cos $&sin (aa-cO + v,, 

(10) 

The equatorial projection of D will be 


rt~D cos di> 

(U) 

so that the fringe rate becomes 


VF=z-(o) c rbVh cos 6 sin (a^-a ) + v 0 

(12) 


For r$=250 km and A =13 cm, the coefficient oi e n- 
/A=140 Hz. Equation (8) yields the time-delay quantity 15 
and is useful when the stations are able to receive wide- 
band signals. For a narrow band signal, as in the present 
invention, Equation (12) is more useful as it yields the 
fringe rate. In otherwords, for a narrow band signal it is 
preferable to operate in the frequency domain, rather 20 
than in the time domain. The frequency offset term v 0 
will not be known a priori with the necessary 50ju,Hz 
(i.e. A///=2X 10- 14 ) accuracy and therefore should be 
simultaneously estimated with the baseline parameters 
tb and cc b. In order for the estimation to be valid, it is 25 
critical that the frequency systems driving the clocks 
have no diurnal dependence because v 0 will then be- 
come confused with the tb and oc b terms. 

From the foregoing, it can be appreciated that the 
fringe rate (fringe frequency observed) as the output of 30 
the cross correlation function may be used to estimate 
the equatorial baseline component Tb and oc b for crustal 
deformation monitoring applications. By measuring r b 
and oc b over an interval of days, months and years, it 
will be possible to detect crustal strain accumulations 35 
and may significantly contribute to an earthquake pre- 
dictive system. 

Although a particular embodiment has been illus- 
trated and described, it is recognized that modifications 
and variations may readily occur to those skilled in the 40 
art such as substitution of radio signals from Global 
Positioning System (GPS). With the advent of the GPS 
as an available signal source several important features 
can be implemented. When fully deployed the GPS will 
consist of 24 satellites in 12 hour period orbits. At least 45 
four satellites will be in view at every global location at 
all times. The satellite orbits will be known from earth 
based tracking data to an accuracy of two meters. Such 
orbital accuracy translates to be 2 cm accuracy in differ- 
ential earth positioning measurements for station sepa- 50 
rations of approximately 120 Km. The GPS radio trans- 
missions are of a broadband type unlike the narrowband 
ALSEP signals discussed previously. By sampling the 
interferometer phase at two or more windows within 
the two bands (1.22 and 1.57 GHz) transmitted by GPS, 55 
it is possible to construct a highly precise delay observ- 
able by the method of bandwidth systhesis described by 
J. B. Thomas, “An Analysis of Long Baseline Radio 
Interferometry,” Part III, Jet Propulsion Laboratory 
Technical Report 32-1526, Vol. XVI, pp. 47-64 (Au- 60 
gust 1973). The two bands separated by 350 MHz pro- 
vides the method for removing the ionospheric error by 
electromagnetic dispersion as descried by R. L. Ko- 
ehler, “Radio Propagation Measurements of Pulsed 
Plasma Streams From the Sun Using Pioneer Space- 65 
craft ” J. Geophys. Res. 73, 15, 4883-4894 (1968). By 
measuring the delay obtained from four or more satel- 
lites at widely separated sky locations, Eq. 8 can be used 


to estimate n, and oc ^ as before, but now the term 
Z^=D sin 8b can be determined giving a full three di- 
mensional baseline description. One of the four delay 
observations is used to estimate the arbitrary phase 
offset term of the bandwidth synthesis, analogous to 
clock synchronization or clock offset term. The near 
realtime advantages of using telephone lines remains as 
in the case of ALSEP because the GPS transmissions 
have high signal strength across their broadcast bands. 

The availability of several satellites simultaneously in 
view begins to approximate the ARIES operational 
method using quasars. The four or more satellite signal 
sources allows the baseline components and clock offset 
to be estimated simultaneously in a time span deter- 
mined by how rapidly the delay observables can be 
derived at the two frequency bands from each satellite. 
Depending on the particular implementation of tele- 
communication elements, the necessary delay measure- 
ments could be collected every few minutes thus signifi- 
cantly reducing the demands upon the frequency and 
timing system. For example, to derive delay precisions 
equivalent to 10 cm every 10 minutes requires a fre- 
quency stability of A///=5x 10~ 13 , readily met by 
commercially available and relatively inexpensive ru- 
bidium frequency systems. Based upon 90 minutes of 
data an effective 3 cm error will be contributed by the 
time and frequency system. 

An important advantage of detecting GPS by cross 
correlation is that no knowledge of the exact nature of 
the transmitted pseudo random digital code is required, 
unlike the case with the original GPS implementation. 
In other words, cross correlation detection of GPS 
transmissions eliminate need for explicit knowledge of 
the form of the digital code being transmitted, which 
may in certain circumstances be transmitted in a secret 
form. Thus, the GPS will become a geodetic resource 
to all users. 

Although particular embodiments of the invention 
have been described and illustrated herein, it is recog- 
nized that modifications and variations may readily 
occur to those skilled in the art. It is therefore intended 
that the claims be interpreted to cover such modifica- 
tions and variations. 

What is claimed is: 

1. In a system for monitoring crustal deformation 
using long-baseline radio interferometry, the method of 
receiving at two independent stations a radio signal of 
limited spectrum from a source at a known location is 
space, maintaining synchronized clocks at said station, 
time tagging the radio signal received at each station, 
transmitting the time-tagged radio signal received at 
each station to a convenient location for cross correla- 
tion using communication channels having a bandwidth 
sufficient for said limited spectrum, adjusting the trans- 
mission delay of one radio signal relative to the other to 
synchronize time tags prior to cross correlation, and 
estimating the interferometer delay or delay rate in 
order to derive baselines for near real-time crustal de- 
formation monitoring. 

2. The method of claim 1 wherein said source is an 
artificial narrow band radio signal, and said communi- 
cation channels are comprised of telephone lines. 

3. The method of claim 2 wherein said source is an 
S-band transmitter on the moon. 

4. The method of claim 2 wherein said source is an 
L-band signal from a global positioning system satellite. 
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5. The method defined by claim 4 including rapid 
simultaneous estimating of baseline using several global 
positioning system satellites in order to significantly 
reduce the stability requirements of the time and fre- 
quency systems. 

6. In a monitoring system, using a pair of antennas at 
separated stations on earth for receiving a radio signal 
of limited spectrum from an artificial source at a known 
location in space and a pair of synchronized clocks, one 
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clock at each station for time tagging the radio signal 
received, the method of transmitting the time tagged 
radio signals from the stations to a third station, adjust- 
ing the transmission delay time of one radio signal rela- 
5 tive to the other to synchronize time tags, and cross 
correlating the signals thus time synchronized for ob- 
taining fringe frequency observations to be used in near 
real-time crustal deformation monitoring. 

* * * * * 
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